Hydrogen in Ag-doped ZnO: theoretical calculations 
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Based on density functional theory calculations, we systematically investigate the behaviors of a H 
atom in Ag-doped ZnO, involving the preference sites, diffusion behaviors, the electronic structures 
and vibrational properties. We find that a H atom can migrate to the doped Ag to form a Ag-H 
complex by overcoming energy barriers of 0.3 - 1.0 eV. The lowest-energy site for H location is 
the bond center of a Ag-0 in the basal plane. Moreover, H can migrate between this site and its 
equivalent sites with energy cost of less than 0.5 eV. In contrast, dissociation of such a Ag-H complex 
needs energy of about 1.1 - 1.3 eV. This implies that the Ag-H complexes can commonly exist in 
the Ag-doped ZnO, which have a negative effect on the desirable p-type carrier concentrations of 
Ag-doped ZnO. In addition, based on the frozen phonon calculation, the vibrational properties of 
ZnO with a Ag-H complex are predicted. Some new vibrational modes associated with the Ag-H 
complex present in the vibrational spectrum of the system. 
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PACS numbers: 66.30.J-, 71.55.Gs, 71.15.Nc 



I. INTRODUCTION 

ZnO is a promising material for short wavelength opto- 
electronic device, due to its various attractive properties, 
such as optical 0, 0, S| , photoelectric [H and piezoelec- 
tric H properties. However, the difficulty in fabricating 
p-type ZnO restricts the application of ZnO in future [6j . 
Therefore, many efforts have been made to achieve p-type 
ZnO, with using many techniques and dopants. Among 
various dopants, either group I elements (Li, Na, and K) 
or group IB elements (Cu, Ag, and Au) are good 

candidates for substitution of Zn. However, with smaller 
ionic radii, group I elements prefer to occupy the intersti- 
tial sites rather than substitutional sites, acting as donors 
(To| . In contrast, group IB elements have suitable size for 
substitution of Zn. So far, the Ag doped p-type ZnO has 
been successfully achieved in the experiment [ll[. 

On the other hand, the role of H in ZnO was paid much 
attention, because H is usually to be unintentionally- 
doped in as-grown ZnO bullc and it always acts as a 
"hole-killer" in ZnO [HI, [H, EH]- Previous studies re- 
vealed that the activation energ y fo r H diffusion in pure 
ZnO was about 0.17 - 0.50 e V [H [H, G3 . which indi- 
cated that a H atom could migrate in ZnO easily at low 
temperature. In ZnO with p-type dopants such as N, 
Cu and Li, H can be easily trapped by the dopants to 
form the impurity complexes. For instance, in N doped 
ZnO, H preferably located at the antibond site of N-Zn 
bond to form a N-H complex [l8|], and the N-H vibra- 
tional modes were observed by Raman spectroscopy in 
experiment [l9[. While in the case of Li doped ZnO, the 
preferred site for H was the bond center of a Li-0 bond 
parallel to the c axis. Such a Li-H complex was studied 
both experimentally and theoretically [Ifl HH H [H, Hi] . 
Except for the Li-H and N-H complexes, the complexes 
of H with transition metals such as Cu, Fe and Ni, were 
also studied [25, 26], in which presence of H significantly 



reduced the effect of the spin polarization. Up to date, 
behavior of H in Ag doped p-type ZnO has not been re- 
ported yet. 

In this paper, based on density functional theory cal- 
culations, we find that a H atom can be trapped by the 
doped Ag atom in ZnO to form a Ag-H complex, and the 
lowest-energy site (ground state site) is the bond center 
of Ag-0 in the basal plane. In addition, the vibrational 
features associated with the Ag-H complex are also dis- 
cussed in this paper. 



II. COMPUTATIONAL DETAILS 

Our calculations are based on the density func- 
tional theory implemented in the SIESTA program [27| . 
Double- £ basis plus polarization (DZP) sets[28( are used 
for all the concerned atoms, and spin polarization is taken 
into account in our calculations. A 72-atom supercell 
consisting of 3 x 3 x 2 primitive unit cells is employed, in 
which a Zn atom is substituted by a Ag atom, and a H 
atom is located at some typical sites around the doped 
Ag respectively. The lattice constants for each case are 
optimized and periodic boundary condition is applied. 
The Brillouin zone is sampled with a set of fc-point grid s 
(3x3x2) according to the Monkhorst-Pack scheme [29l ]. 
The configurations are fully relaxed with using the con- 
jugate gradient method until the Hcllman-Fcynman force 
on each atom is less than 0.02 eV/A. 

With these settings, we optimize the geometry of the 
perfect wurtzite ZnO (w-ZnO) both with the local den- 
sity approximation (LDA) and with the Perdew-Burkc- 
Ernzerhos generalized gradient approximation (GGA- 
PBE) . The obtained optimal crystallographic parameters 
of w-ZnO from LDA calculations are a— 3.278 A, c=5.237 
A, and w=0.382, which are slightly larger than those ( a 
=3.25 A, c=5.207 A, and u=0.382) measured from the 
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experiment [30] at room temperature. It is noted that the 
previous experiment demonstrated that the lattice con- 
stants of ZnO decreased slightly when the temperature of 
the system increased [3lJ. Therefore, our calculated pa- 
rameters are consistent with the experimental values at 
low temperature. In contrast, the lattice constants eval- 
uated at the level of GGA are larger than those at the 
level of LDA. On the other hand, the formation enthalpy 
(-3.43 eV/ZnO) of w-ZnO obtained from the LDA calcu- 
lation does also match that (-3.61 eV/ZnO) derived from 
the experiments. Therefore, we perform our calculations 
at the level of LDA. In addition, we test the convergence 
of the super cell size by calculating the defect formation 
energy of the Agzn in the supercells of 2x2x2, 3x2x2 
and 3x3x3 respectively, in which the lattice constants 
are optimized in each case. Here the defect formation 
energy is defined as the difference of the total energies 
between the defect case and the corresponding perfect 
case. We find that the difference of the defect formation 
energies between 2x2x2 and 3x2x2 supercells is as 
large as 0.5 eV, whereas for the supercells of 3 x 3 x 2 and 
3x3x3, the difference of the defect formation energies 
is about 0.02 eV only. This shows that the interaction of 
the doped Ag with its images is very weak in supercell 
3x3x2, and thus the influence of the images of Ag can 
be neglected. In addition, the obtained bond lengths as 
well as the electronic structures in 3 x 3 x 2 are almost 
the same as that in 3 x 3 x 2 supercell. So the results 
achieved from the 3x3x2 supercell are reliable. 

The frozen phonon approximation [32| is employed to 
explore the vibrational properties of the typical systems 
with Ag-H complexes. In our calculations, the atoms 
in each of the systems are displaced one by one from 
their equilibrium positions along three Cartesian direc- 
tions and the reverse directions, with an amplitude of 
0.04 Bohr. By the numerical derivatives for the displace- 
ments of each atom, the force constants are obtained, 
which are used to build up the dynamical matrix of this 
system. By solving the dynamical equation 

UI 2 MiUi, a = ^ Ci,a;j,pUj,p, (1) 

the vibrational frequencies uj and the corresponding 
eigenmodes yielded, where the Mj is the mass of 

ith atom, Ci iCt; j t p is the force constant, and a (/3) means 
the direction of x, y or z. 



III. DIFFUSION OF A H ATOM IN AG-DOPED 

ZNO 

Firstly, we replace a Zn atom in the perfect ZnO su- 
percell by a Ag atom, then optimize the system with the 
considerations above. Then a H atom is located around 
some typical O atoms, which are notated as Oi, On, 
Oiu, Oiv and Oy respectively, as shown in Fig.l. These 



TABLE I: The optimized coordinates of H at the stable 
(metastable) sites near the doped Ag. The coordinates of the 
doped Ag is also listed. The positions for H refer to Fig.l. 
The lattice vectors of the supercell are a=(9.85, 0.0, 0.0), 
6=(-4.93, 8.53, 0.0), and c=(0.0, 0.0, 10.65). The unit of the 
length is in angstrom. 





Al 


A2 


A3 


Bl 


B2 


B3 


B4 


B5 


B6 


B7 


Ag 


X 


1.71 


1.64 


-1.13 


1.65 


0.63 


-0.58 


1.07 


-0.03 


2.68 


1.03 


1.64 


y 


2.87 


2.84 


4.65 


2.86 


3.43 


4.13 


3.14 


5.55 


3.42 


5.96 


2.85 


z 


8.80 


3.37 


4.18 


6.62 


4.91 


3.86 


1.92 


6.62 


4.92 


4.57 


5.32 



O atoms are near the doped Ag within the fifth neigh- 
bour. For the structure of ZnO bulk, there are two kinds 
of inequivalent Zn — O bonds which is either parallel to c 
axis or within the basal plane respectively. Therefore, for 
the Ag-doped ZnO bulk, there are four inequivalent sites 
available for the location of H around each concerned O 
atom: bond center (BC) sites and antibond (AB) sites 
of Zn(Ag) — O bonds parallel to the c axis and those 
of Zn(Ag) — O bonds in the basal plane. From the cal- 
culations, we find that some of these sites are unstable 
for H, where the located H atom spontaneously moves to 
other places during full relaxation. For example, when a 
H atom is located at either AB site of Zn — On bond 
parallel to the c axis or BC site of Zn — On bond in 
the basal plane, it spontaneously migrates to a BC site 
or AB site of Ag — On bond during geometry optimiza- 
tion. For convenience, some stable or metastable sites 
for H are marked with characters of A n (n=f to 3) for 
AB sites and B n (n=l to 7) for BC sites in Fig.f . The 
coordinates of these sites and the doped Ag are listed in 
Table I. 

Among these sites, site B2 in Fig.l is the ground state 
site for H location, which is similar to that of H in Cu 
doped ZnO [33]. Meanwhile, site ^43 is another preference 
site for H, where the energy of the system is only 0.03 eV 
higher than that at site B2. Nevertheless, the obtained 
energies for the other metastable sites are higher than 
that at site B2 by less than 1.0 eV. 

Since the ground state and metastable sites for H are 
obtained, we now study the diffusion behaviors of H near 
the doped Ag. Eleven diffusion paths are considered, 
which are marked in Fig.l. By using the climbing image 
nudged elastic band (CI-NEB) scheme [H], [Uj], the ener- 
gies for the images in each diffusion path are calculated. 
Taking the lowest energy of the system (with H at site 
B2) as a reference, the energy profiles for each concerned 
diffusion path are plotted in Fig. 2. As shown in Fig. 2 
(a), once a H atom arrives at the site of either Al or Bl 
near Oj, it can migrate to the ground state site (B2), by 
overcoming energy barriers of less than f .0 eV. 

In contrast, when a H atom reaches site B4 near Om, 
it can diffuse to site B2 via site B3 (path-4, path-3), with 
energy cost of less than 0.5 eV. For a H atom at site A2, it 
prefers to migrate towards site B2 directly through path- 
5, rather than towards site AZ via site B3 (path-6, path- 
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FIG. 1: (Color online) Schematic showing of the concerned 
diffusion paths of H around the doped Ag in ZnO, indicated 
by path-1 to path-11. As labeled, green, purple, red and gray 
balls represent O, Zn, Ag and H atoms respectively. Some O 
atoms near the Ag atom are marked with Romanic numbers. 
The stable and metastable sites for H location in the diffusion 
paths are marked with symbols of A n (n=l to 3) and B n (n=l 
to 7). 

7), because the energy barrier in the former case is about 
0.4 eV, while that in the latter case is about 0.7 eV, which 
can be seen in Fig. 2 (b). Meanwhile, if a H atom reaches 
near Ov (at site Bh) or Oiy (site B7), it can reach site 
A3 through path-8 or path-11, with overcoming a low 
barrier of about 0.4 or 0.3 eV, as shown in Fig. 2 (b,c). 
Interestingly, the energy barriers for H diffusing between 
the equivalent sites of B2 and B6 are less than 0.5 eV 
(Fig. 2 (c), path-9), which indicates that a H atom can 
move "freely" between B2 and its equivalent sites near 
the doped Ag at room temperature. Such local motion of 
H near the doped Ag was also found in the Cu-doped ZnO 
[111 . Moreover, H can diffuse between the low energy sites 
of B2 and ^43 through two pathways. One is from site B2 
to site A3 directly via path-10 with energy barriers of less 
than 0.9 eV, and the other is to site B3 through path- 
3 firstly, then to site A3 through path-7 by overcoming 
the energy barriers of no more than 0.7 eV. On the other 
hand, one can find from the Fig. 2 that once a H atom 
locates at either site B2 or A3, the dissociation of such 
a Ag-H complex needs the energy cost as high as about 
1.1- 1.3 eV. 

From above, we can conclude that a H atom can diffuse 
to the doped Ag to form a Ag-H complex in ZnO through 
many paths, and the lowest diffusion barrier is 0.3 eV 
only. In contrast, the dissociation of a Ag-H complex 
requires energy costs of more than 1.1 eV. This indicates 
that the Ag-H complex may easily exist in the Ag-doped 
ZnO. We stress that such stable complexes of Ag-H in 
ZnO play a negative role in p-type ZnO with doped Ag. 

To go further, we pay our attention to the influence 
of a Ag-H complex on electronic structures of the sys- 
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FIG. 2: Energy profiles calculated with CI-NEB for the diffu- 
sion paths of H in Ag-doped ZnO: (a) from path-1 to path-4, 
(b) from path-5 to path-8, and (c) from path-9 to path-11. 
The calculated relative energies for the images in each path 
are indicated with stars. The diffusion paths and reaction 
coordinates shown in this figure refer to Fig.l. 



tern. The calculated total density of states (TDOS) for 
the systems with a H atom at site B2 and A3 are shown 
in Fig. 3 a and b respectively. For comparison, the TDOS 
for the Ag monodoped ZnO and the perfect ZnO with 
the same size are shown in Fig. 3 (c). Obviously shown in 
Fig. 3 (a and b), some occupied defect states appear near 
Fermi level for the cases of H at site B2 and site A3. By 
analysis of the local density of states (LDOS), we find 
that these defect states mainly come from the contribu- 
tion of the doped Ag atom and its neighbouring O atoms, 
and little from H and its bonded O atom, as displayed 
in Fig. 3 (d and e). Furthermore, the projected density 
of states (PDOS) analysis reveals that these states are 
mainly contributed from 3d orbitals of the Ag atom and 
2p orbitals of the O atoms. 
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FIG. 3: (Color online) The total density of states (TDOS) of 
the system for (a) the Ag-H complex with H at site B2, (b) the 
Ag-H complex with H at site A3, (c) ZnO with and without Ag 
doping, and the corresponding local density of states (LDOS) 
(d), (e) and (f). The dash line in (c) stands for the TDOS of 
the perfect ZnO, and the solid lines in (c) stand for Ag-doped 
ZnO with majority spin and minority spin. Ol refers to the 
O atom being far away from the the doped Ag; 02 and 03 
stand for the O atoms in Ag-O parallel and perpendicular to 
the c axis respectively, and 04 for the O atom in H-O bond. 
The dot-dash lines indicate the Fermi levels, which are shifted 
to be zero. 



IV. VIBRATIONAL PROPERTIES 

Vibration is a fundamental property of a system, which 
is closely coupled with the structure of the system. For 
such a Ag-H complex in ZnO, some new vibrational fea- 
tures should appear in the vibrational spectrum. With 
the frozen phonon approximation, we explore the vibra- 
tional properties of the ZnO containing a Ag-H complex. 
As a reference, we firstly calculate the total vibrational 
density of states (TVDOS) for the perfect ZnO, which is 
plotted in Fig. 4 (b). The high frequency range for the 
perfect ZnO is below 570 cm -1 , and a silent region exists 
between 270 and 400 cm -1 . This is consistent with the 
previous report [36| . 

For the system containing a Ag-H complex with H at 
site B2, a typical stretching mode of H-0 bond is found 
to be at 3275 cm -1 , and no state between 782 and 3275 
cm,- 1 . The calculated TVDOS in the region below 782 
arT 1 is shown in Fig. 4 (b). From Fig. 4 (b), we find that 
some new vibrational states, which are marked with A, 
B, C and D respectively, appear in the silent frequency 
region of the perfect ZnO. 

The local vibrational density of states (LVDOS) anal- 
ysis reveals that the states at lower frequency region of 
TVDOS (below 270 cm,- 1 ) for the perfect ZnO comes 
from the contribution of the Zn atoms, and those at 
higher frequency region (above 400 cm -1 ) from the O 
atoms. This remains in the TVDOS for ZnO containing a 
Ag-H complex. Furthermore, we find that the vibrational 



FIG. 4: (Color online) (a) The local vibrational density of 
states (LVDOS) and (b) total vibrational density of states 
(TVDOS) for the Ag-H complex with H at site B2, where 
frequency range is below 800 cm - . The stretch modes of H- 
O are not included. 1,2,3,4 and 5 in (a) represent the averaged 
LVDOS for the O atoms far from Ag, the O atoms around Ag, 
the O atom bonding with H, Ag and H respectively. In (b), 
the red dash line stands for the TVDOS of the perfect ZnO, 
and the solid line for that of the complex system with Ag-H. 
The states marked with the characters indicate the new states 
induced by the Ag-H complex. 



states marked with characters of A (below 400 cm -1 ) and 
B (above 570 cm -1 ) in Fig. 4 (b) are mainly from the con- 
tribution of four nearest neighbouring O atoms around 
Ag, and little from the O atoms far from Ag, which can be 
seen obviously in Fig. 4 (a) . The presence of these modes 
are ascribed to the local structural distortion induced by 
the Ag-H complex. By examining the configuration of the 
system, we find that the Ag-0 bonds both parallel and 
perpendicular to c axis are all enlarged by about 0.2 A. 
Consequently, the Zn-0 bonds in the vicinity are short- 
ened or enlarged within 0.04 A, and the bond lengths of 
the next neighbouring Zn-0 are altered slightly. Such 
structural changes result in presence of these new states 
near A and B in Fig. 4 (b). Meanwhile, the vibrational 
states marked with C (at about 646 cm- 1 ) and D (782 
cm^ 1 ) are attributed to the bending modes of the H-0 
bond. In addition, the states contributed from the doped 
Ag mainly emerge at the low frequency region, as shown 
in Fig. 4 (a). 

For the case of a H atom at site ^43, we also calculate 
the TVDOS of the system, and its main features are sim- 
ilar to those of H at site B2. The stretching mode of the 
H-0 bond locates at the frequency of 3287 cm -1 , and the 
bending modes of the H-0 bond are at 768 and 879 cm" 1 
respectively. Similarly, some new states associated with 
the neighbouring O atoms around the doped Ag emerge 
near the frequency regions of below 400 and above 570 
cm- 1 respectively. Especially, a mode corresponding to 
displacement of H and its bonded O atom in the same 
direction is found at the frequency of 273 cm- 1 . These 
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typical vibrational modes may drip a hint to detect the 
presence of such a Ag-H complex in experiment. 

V. SUMMARY 

In summary, the site preferences and diffusion behav- 
iors of a H atom in Ag-doped ZnO are investigated based 
on density functional theory calculations. We find that 
the lowest-energy site for H in the Ag-doped ZnO is bond- 
centered sites of Ag-0 in the basal plane. H can migrate 
to these sites through some typical diffusion paths with 
energy barriers of 0.3 - 1.0 eV, most of which are less 
than 0.7 eV. H can diffuse between its ground state sites 
with the energy cost of less than 0.5 eV. In contrast, 
releasing of this H atom from the doped Ag requires 
energy of about 1.1 - 1.3 eV, which indicates that the 
Ag-H complex may commonly exist in Ag-doped ZnO. 



By calculating the vibrational properties of ZnO with a 
Ag-H complex, we find some new vibrational modes in 
the silent frequency region of perfect ZnO, which are at- 
tributed to the distortion of the O atoms nearby induced 
by the Ag-H complex. 



VI. ACKNOWLEDGMENTS 

This work is supported by the University of Science 
and Technology of China, the Chinese academy of Sci- 
ence, National Science Foundation of China (Grant Nos. 
NSFC10574115 and NSFC50721091). B. C. Pan thanks 
the support of National Basic Research Program of China 
(2006CB922000). This indicates that the Ag-H com- 
plexes can commonly exist in ZnO. The HP-LHPC of 
USTC is acknowledged for computational support. 



[1] D. C. Look, D. C. Reynolds, J. R. Sizelove, R. L. Jones, 

C. W. Litton, G. Cantwell, and W. C. Harsch, Solid State 

Commun. 105, 399 (1998). 
[2] S. Y. Lee, Y. Li, J.S. Lee, J.K. Lee, M. Nastasi, S.A. 

Crooker, Q.X. Jia, H.S. Lee, and J.S. Kang, Appl. Phys. 

Lett. 85, 218 (2003). 
[3] H. S. Kang, J. S. Kang, J. W. Kim, and S. Y. Lee, J. 

Appl. Phys. 95, 1246 (2004). 
[4] W. W. Wenas, A. Yamada, and K. Takahashi, J. Appl. 

Phys. 70, 7119 (1991). 
[5] T. Mitsuyu, S. Ono, and K. Wasa, J. Appl. Phys. 51, 

2646 (1980). 

[6] S. B. Zhang, S. H. Wei, and A. Zunger, J. Appl. Phys. 

83, 3192 (1998), and references therein. 
[7] A. Zunger, Appl. Phys. Lett. 83, 57 (2003). 
[8] S.B. Zhang, S.-H. Wei, and A. Zunger, Phys. Rev. B 63, 

075205 (2001). 
[9] S.-H. Wei, Comput. Mater. Sci. 30, 337 (2004). 
[10] D.C. Look, R.L. Jones, JR. Sizelove, NY. Garces, N.C. 

Giles, and L. E. Halliburton, Phys. Status Solid A 195, 

171 (2003). 

[11] H.S. Kang, B.D. Ahn, J.H. Kim, G.H. Kim, S.H. Lim, 
H.W. Chang, and SY. Lee Appl. Phys. Lett. 88, 202108 
(2006). 

[12] C. G. Van de Walle, Phys. Rev. Lett. 85, 1012 (2000). 

[13] S. F. J. Cox, E. A. Davis, S. P. Cottrell, P. J. C. King, J. 
S. Lord, J. M. Gil, H. V. Alberto, R. C. Vilao, J. Piroto 
Duarte, N. Ayres de Campos, A. Weidinger, R. L. Lichti, 
and S. J. C. Irvine, Phys. Rev. Lett. 86, 2601 (2001). 

[14] D. M. Hofmann, A. Hofstaetter, F. Leiter, H. Zhou, F. 
Henecker, B. K. Meyer, S. B. Orlinskii, J. Schmidt and 
P. G. Baranov, Phys. Rev. Lett. 88, 045504 (2002). 

[15] K. Ip, M. E. Overberg, Y. M. Heo, D. P. Norton, S. J. 
Pearton, C. E. Stutz, B. Luo, F. Ren, D.C. Look, and 
J.M. Zavada, Appl. Phys. Lett. 82, 385 (2003). 

[16] N. H. Nickel, Phys. Rev. B 73, 195204 (2006). 

[17] M.G. Wardle, J. P. Goss, and PR. Briddon, Phys. Rev. 
Lett. 96, 205504 (2006). 

[18] X. N. Li, B. Keyes, S. Asher, S. B. Zhang, S. H. Wei, T. 
J. Coutts, S. Limpijumnong, C. G. Van de Walle, Appl. 



Phys. Lett. 86, 122107 (2005). 
[19] O. S. Kumar, E. Watanabe, R. Nakai, N. Nishimoto, and 

Y. Fujita, J. Cryst. Growth 298, 491 (2007). 
[20] L. E. Hallibuton, L. J. Wang, L. H. Bai, N. Y. Garces, 

N. C. Giles, M. J. Callahan, and B. G. Wang, J. Appl. 

Phys. 96, 7168 (2004). 
[21] M. G. Wardle, J. P. Goss, and P. R. Briddon, Phys. Rev. 

B 71, 155205 (2005). 
[22] G. A. Shi, M. Stavola, and W. B. Fowler, Phys. Rev. B 

73, 081201(R) (2006). 
[23] K. R. Martin, P. Baney, G. Shi, M. Stavola, and W. B. 

Fowler, Phys. Rev. B 73, 235209 (2006). 
[24] E. V. Lavrov, F. Borrnert, and J. Weber, Phys. Rev. B 

71, 035205 (2005). 
[25] M. G. Wardle, J. P. Goss, and P. R. Briddon, Phys. Rev. 

B 72, 155108 (2005). 
[26] F. Borrnert, E. V. Lavrov, and J. Weber, Phys. Rev. B 

75, 205202 (2007). 
[27] D. Sanchez-Portal, P. Ordejon, E. Artacho and J. M. 

Soler, Int. J. Quantum Che 65, 453 (1997); N.Troullier 

and J. L. Martins, Phys. Rev. B 43, 1993 (1991). 
[28] J. M. Soler, E. Artacho, J. D Gale, A. Garcia, J. Jun- 

quera, P. Ordejon, and D. Sanchez-Portal, Phys. Con- 
dens. Matter 14, 2745 (2002), and references therein. 
[29] H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 

(1976). 

[30] J. Albertsson, S.C. Abrahams, and A. Kvick, Acta Crys- 
tallogr., sect. B: Struct. Sci. 45, 34 (1989); Handbook of 
Chemistry and Physics, 85th ed., edited by DR. Lide 
(CRC Boca Raton, FL, 2004). 

[31] L. Yan, C.K. Ong, X.S. Rao, J. Appl. Phys. 96, 508 
(2004). 

[32] M.T. Yin, ML. Cohen, Phys. Rev. B 26, 3259 (1982). 
[33] J. Hu, and B.C. Pan, J. Phys. Chem. C 112, 19142 (2008). 
[34] G. Henkelman, B. P. Uberuaga, and H. Jonsson, J. Chem. 

Phys. 113, 9901 (2000). 
[35] G. Henkelman, and H. Jonsson, J. Chem. Phys. 113, 

9978 (2000). 

[36] A. W. Hewot, Solid State Commun. 8, 187 (1970). 



